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Abstract. We exhibit seven linear codes exceeding the current best known minimum distance 
d for their dimension k and block length n. Each code is defined over Wg, and their invari- 
ants [n,k,d] are given by [49,13,27], [49,14,26], [49,16,24], [49,17,23], [49,19,21], [49,25,16] 
and [49,26,15]. Our method includes an exhaustive search of all monomial evaluation codes 
generated by points in the [0, 5] x [0, 5] lattice square. 



1. Introduction 

The basic invariants of a linear code C over the finite field ¥ q are its dimension k and its block 
length n. The code is the image in IF" of a k x n matrix M over ¥ q . The minimum distance 
d of C is the smallest Hamming weight of any nonzero linear combination of the rows of M. 
These invariants are conventionally recorded as a triple [n,k,d]. One usually seeks codes with 
minimum distance d as large as possible for a given block length n and dimension k. There are 
theoretical upper bounds for the minimum distance, although in many cases the largest known 
example falls short of these bounds. Grassl |Graj catalogues the best known linear codes in this 
sense, and we refer to any code with larger minimum distance as a champion code\j\. 

In this paper we find seven new champion codes defined over ¥$ by considering the class of 
generalised toric codes introduced by Little [Litll] (see $2]). 

Theorem 1.1. There are precisely five [n, k,d] profiles of champion generalised toric codes over 
¥g generated by collections of points in a [0,5] x [0,5] square. There are (at least) another two 
champion codes arising from a [0,6] x [0,6] square. These profiles are listed in Table [ZJ 

The codes we consider correspond to subsets of lattice points in the [0, m] x [0, m] square for 
small values of m. In principle, our approach is to enumerate all such sets of points, construct 
the corresponding code over various small fields and then compute its invariants. In $3] we 
describe an algorithm to enumerate all such sets of points up to affine equivalence. To prove the 
theorem, we implement and run this algorithm to completion for m = 5, and a partial search 
when m = 6. Since our algorithm is exhaustive, we can say more: 

Theorem 1.2. There are no champion generalised toric codes over ¥ q when q < 7. 

It is worth noting that the new Bound A of [GBS12J docs not improve the existing theoretical 
upper bounds for minimum distance in the range of the champion codes we discover, so we 
cannot say that the codes of Theorem 11.11 are the best possible amongst all linear codes. 
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^The profile [36, 19, 12] over F7 discovered in [BK12] is not currently recorded in |Gra| . but we will not regard 
any code matching these invariants as a champion. 
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Table 1. The [n,k, d] invariants of new champion generalised toric codes over 
Fg. In each case a single example is given, contained in the smallest possible 
[0, m] x [0,m] square, and illustrated in Figure [TJ 



2. Generalised toric codes 

Recall from [HanOOj that a toric code C, over a sufficiently large field ¥ q , is determined by a 
convex lattice polygon P C 1? <g> Q as follows. Suppose that P lies in a [0, m] x [0, m] square, 
where q > m + 2. Then C is given by the image of a k x n matrix M whose rows are generated 
by evaluating each lattice point (a, b) £ P n Z 2 (regarded as a monomial x a y h ) at each vector 
of the torus (F*) 2 . The dimension of C is the number of rows k = \P n 1? | of M and the block 
length is the number of columns n = (q — l) 2 . 

A generalised toric code is constructed in the same way as a toric code but with the possibility 
of omitting one or more of the lattice points of P; equivalently, one may remove rows from the 
generating matrix M. Conceptually this allows one to delete any particularly short vectors in 
C that arise as rows of M. Although of course this does not necessarily increase the minimum 
distance, they have recently also been proving fruitful in the search for champions |Rua091 
IAHV091 ICMO091 ILitllj . Little introduced the study of generalised toric codes, and found a 
champion [49, 12, 28] code over Fs coming from a particular subset of points of a polygon in a 
[0,5] x [0,5] square [Litll] : this is illustrated in Figure EJa). Systematic attempts to produce 
examples of champion generalised toric codes have been performed over F4 and F5 |AHV09) . 
and over F 7 , F 8 and F 9 |CMU09| . 
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Figure 1. Example point configurations for the new champion generalised codes 
over Fg listed in Table [TJ 



Small polygons. In |BK12j we assembled a comprehensive database of lattice polygons that 
are contained, up to lattice automorphism and translation, in a [0, m] x [0, m] square, for m < 7. 
This database can be interrogated online via the Graded Ring Database [BKaj, or from within 
the computational algebra software Magma [BCP97| . In [BK12] we checked the toric code 
corresponding to each such polygon, over all prime-powered fields ¥ q for m + 2 < q < 9, and 
found a single new champion code, defined over F7 with invariants [36, 19, 12]. In this paper we 
consider all generalised toric codes associated with polygons that lie in a [0, m] x [0, m] square, 
for m < 5 and m + 2 < q < 9. We also present partial results when m = 6; in this case 
the number of possible generalised toric codes is too large to be searched systematically using 
current techniques. The invariants listed in Table [T] are those of the champion codes we found 
with this search. 




FIGURE 2. (a) The [49, 12,28] code over F 8 described by Little (Litllj : (b) One 
of 448 non-equivalent ways of generating a [49, 12, 28] code over Fs from points 
in a [0, 4] x [0, 4] square. 



Multiplicity of champions. Champion profiles such as [49, 13, 27] over Fs are often achieved 
in many non-equivalent ways: this case, for example, is realised by four different sets of points in 
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a [0,4] x [0,4] square, distinct up to lattice automorphism and translation; see Figure [3l Little's 
champion [49, 12, 28] code is more striking still: we have 448 non-equivalent generalised codes 
in a [0,4] x [0,4] square that yield these invariants. An example is illustrated in Figure [2jb). 



Figure 3. The possible choices of points in a [0,4] x [0,4] square giving a 
[49,13,27] code over F 8 . 

Toric codes often achieve best known minimum distance. We extracted the [n, k, d] 
profiles of all generalised toric codes that matched, or exceeded, the largest minimum distances 
available in Grassl's catalogue over ¥ q for q < 7. All such codes are contained in a [0, 5] x [0, 5] 
square, and so our results are complete. The bounds attained are recorded in Appendix |Aj 
generalised toric codes achieve (or exceed) the current best known minimum distance in 28 of 
the 57 cases. 



The proof of Theorems 11.11 and 11.21 is by a systematic computer search. The hurdle to 
overcome is the sheer number of codes and the complexity of computing their minimum distances 
(often these calculations would take millions of years if the code achieved its apparent minimum 
distance). We describe an algorithm that carries out the enumeration of subsets of lattice points 
of a lattice square up to affine equivalence. We have implemented this in Magma exactly as 
described here, making use of the convex polytopes package [BBKllj . Our code is available to 
download from |BKbj . 

Let S C 7? be a collection of points contained in a [0, m] x [0, m] square generating the 
generalised toric code Cs over ¥ q . For any lattice translation S' = S — u of the points S, 
or more generally for any affine linear transformation s i— >• (s — u)M of S, where u € 71? 
and M € GL2(Z), the toric codes Cs and Cs' are monomially equivalent |LS07| Theorem 4]. 
Therefore it is enough to consider points S up to lattice translation and change of basis. Consider 
now the lattice polygon P := conv(S') C 77? (8> Q. Up to equivalence, P can be assumed to be 
one of the polygons constructed in |BK12j . This motivates our approach: For each polygon P 
contained in a [0, m] x [0, m] square we will generate, up to equivalence, all possible subsets S 
of points P fl 7? such that S contains the vertices vert(P) of P. Insisting that vert(P) C S is a 
natural restriction; if this were not so, Q := conv(S') is a lattice polygon distinct from P in the 
same [0, m] x [0, m] square, and, up to equivalence, that polygon will be considered separately. 

Step 1: Compute the affine automorphism group of P. The first step is to compute the 
affine automorphism group G := AffAut(P) of the polygon P, where any element g £ G can be 
written as a combination of elements of G\^2{7L) and translations. Embed P at height 1 in the 




3. The algorithm 
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lattice Z 3 = Z 2 x Z; for example, simply append a coordinate 1 to each vertex of P. We refer 
to this embedded image of P also as P. Form the cone a with vertex the origin generated by 
the points of the embedded image of P. Let Aut(cr) be the linear automorphism group of a. 
The action of this group on a restricts to a faithful action on P, realising the full group of affine 
lattice automorphisms of P. 

Step 2: Extend the action to subsets of points of P. Clearly v ■ G C vert(P) for any 
v G vert(P). Let V := P n Z 2 \ vert(P) be the set of non- vertex lattice points in P. We fix an 
order of the points of V = {v±, . . . ,Vk}, where k := \V\; G acts on V via permutation. Choose 
a largest orbit 0\ of this action on V . Without loss of generality we assume that 0\ permutes 
the first h\ < k elements. 

Step 3: Enumerate subsets of a largest orbit up to the action. From now onwards we 
regard G as acting on {0, l} k , where an element s G {0, l} k corresponds to a choice of points 
SCPflZ 2 via the obvious map 

s = (&i, . . . , b k ) ^ S = vert(P) U {vi G V | 1 < % < k, h = 1}. 

Let Si be the set of all sequences in {0, up to this action. Enumerating S\ for small 
polygons is straightforward: in particular, S± contains a unique sequence of zeros and also a 
unique sequence containing a single 1. Larger numbers h > 2 of nonzero coefficients depend on 
the /i-transitivity of the action: the number is determined by the orbit-stabiliser theorem, but 
even without applying that the numbers are small enough simply to run through all possibilities 
rejecting those already seen (up to the action). We only compute the results for at most |_A;i / 2j 
nonzero coefficients, since the remaining possibilities are obtained via symmetry by exchanging 
and 1. 

Step 4: Extend subsets to all of V. For each s G S\ we extend s to an element of {0, l} k 
as follows. Let G s := {g G G \ s ■ g = s} < G, and let V s := G V \ k\ < i < k}. By induction 
we can construct the set of possible {0, l} fc ~ fcl -sequences S s , up to the action of G s . Let sU s' 
denote the concatenation of s G {0, l} fcl with an element s' G {0, so that sUs'e {0, l} k . 

Then 

S = {sUs' | s G Si, a' G S s } 
corresponds to the set of all points in P n Z 2 \ vert(P), up to the action of G. 

Example 3.1. Consider the polygon B equal to the [0, 2] x [0, 2] square. This contains 3 2 points, 
and hence there exist 2 32 subsets of points in B. If one insists that vert(-B) C 5 for any subset 
S of points, then this is reduced to 2 3 ~ 4 = 32 possibilities. Considering subsets S only up to 
equivalence, we obtain the twelve equivalence classes illustrated in Figure HI 

At this point, we have a list of all possible subsets of points of P up to affine automorphisms 
(and including its vertices). It remains to identify any champion generalised toric codes. 
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Figure 4. The twelve equivalence classes of points S in the [0,2] x [0,2] square 
B such that vert (5) C 5. 



Step 5: Run a trial minimum distance algorithm to exclude most cases. Fix a prime 
power q and an integer < m < q — 1. Given a set of points S C [0,m] x [0, m], compute the 
generalised toric code C = Cg(F g ). (In our case, S* is one of the subsets of points of a polygon 
enumerated in Step 4.) Let M = Mq be the generator matrix of C. Compute the Hamming 
length of the shortest nonzero vector that is a linear combination over ¥ q of up to four rows of 
M - this is q A — 1 vectors, which is manageable for the small values of q we consider. If this 
length is strictly greater than the best known minimum distance for codes of dimension \S\ and 
block length q — 1, then keep the pair S, q for consideration at the next step; otherwise discard 
the pair S, q, since the corresponding code cannot be a new champion. 

Step 6: Compute the minimum distance of successful trial candidates. In principle 
this is now the hard part. We employ Magma's MinimumDistance function; at heart, it must 
simply check the Hamming weight of all words in the code, and so can easily be expected to 
take too long for our purposes. In practice when m < 5, for the hundreds of thousands of codes 
that pass our four-line trial in Step 5, we can always calculate the minimum distance within 
a few hours, and there are several hundred cases realising champion invariants. So for m < 5 
there are no numerical profiles [n, d, k] other than the five listed in Table Q] that admit champion 
generalised codes. 

Remark 3.2. To get an idea of the effect of symmetries, Table [2] lists the number of subsets of 
points 2( m+1 ) 2 ~ 4 of the [0 , m] x [0, m] square (excluding the vertices, since we handled vertices 
of polygons separately) and the number of such subsets up to equivalence. Of course there are 
eight symmetries of the square, and 8 x 265488 = 2123904; the full group of symmetries of P 
has been exploited. 
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Table 2. A comparison of the number of subsets of points in a [0, m] x [0, m] 
square, and the number of subsets up to equivalence. 
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It is conceivable that there are other 'hidden' symmetries that preserve the invariants of the 
toric code, even if they do not preserve the lattice subset, and the multiplicity of champions 
may be hinting at this. 

Acknowledgments. Our thanks to Tom Coates for several helpful ideas, to John Cannon for 
providing Magma for use on the Imperial College mathematics cluster, and to Andy Thomas 
for technical assistance. The second author is supported by EPSRC grant EP/I008128/1. 
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Appendix A. Generalised toric codes over ¥ q , q < 7 

The largest minimum distance dt achieved by a generalised toric code with given block length n 
and dimension k over the field ¥ q are documented in the following table. For comparison, we also 
give the current best minimum distance d g amongst all linear codes as given by Grassl [Gra]. The 
penultimate column records whether dt > d g , and the final column gives an example collection of 
points realising these invariants. The only example for which a generalised toric code is required 
is [36,7,23] over F 7 . 
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